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ABSTRACT 

The white dwarf (i.e., the primary) in the AE Aqr binary star system is observed 
to spin down at a steady time rate Pnow = 5.64 x 10~^^ss~^, while at the same time 
its UV (£(7y) and X-ray accretion luminosities seem to remain almost unchanged. 
This is a contradiction, however, since the classically estimated spin-down power due 
to decrease of rotational kinetic energy T, T ~ — 10^^ ergs" ^, exceeds the accretion 
luminosity of the primary by a factor ^ 10^ and, as a dominating power, should lead 
either to observable luminosity changes or to other detectable effects. This so-called 
"spin-down problem" can be relaxed under the assumption that the primary is now in a 
phase of rapidly decreasing its differential rotation under constant angular momentum, 
undergoing a nonaxisymmetric transition — i.e., turning over its magnetic symmetry 
axis with respect to its angular momentum axis and eventually becoming a perpendic- 
ular rotator — from a "differential rotation state" (DRS) to a "rigid rotation state" 
(RRS) on a relatively short "nonaxisymmetric DRS-to-RRS transition timescale". If 
so, then the estimated spin-down power, — T ^ 4 x lO^^ergs"^, becomes comparable 
to the observed luminosities and the spin-down problem is drastically simplified. In 
this paper, we present a detailed numerical study of such a "turn-over scenario" , which 
study mainly points out the fact that an observed large spin-down time rate does not 
always imply a large spin-down power. 



Subject headings: stars: individual (AE Aqr) — stars: magnetic fields — stars: rotation 
and spin-down — stars: white dwarfs 



Introduction 



AE Aqr has been studied extensively in the last years (see, e.g.. 


Welsh, Home, & Gomer 


1993 


; Welsh, Home, & Oke 1993|; de Jager et al. 1994, hereafter dJ&94; 


Eracleous et al. 1994, 


hereafter E&94; Meintjes et al. 1994 


, hereafter M&94; Patterson 1994, 


hereafter P94; 


[khsanov 



1995, hereafter 195; Welsh, Home, & Gomer 1995, hereafter WHG95; Casares et al. 1996, hereafter 
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C&96; 


Eracleous & Horne 1996, hereafter EH96; 


[khsanov 1997 


hereafter 197; 


Kuijpers et al. 


1997 


, hereafter K&97; 


Wynn, King, & Horne 1997 


, hereafter WKH97; 


[khsanov 1998 


, hereafter 


198; Welsh, Horne, & Corner 1998|; |Choi, Dotani, & Agrawal 1995 ; 


Ikhsanov 1999, hereafter 199). 



AE Aqr has been classified (see, e.g., P94, § 2.4.2) as member of the DQ Her or Intermediate 
Polar (IP) cataclysmic variables (CVs). Such CVs are binary star systems in which a rotating 
magnetic white dwarf (playing the role of "primary" ) and a late-type main-sequence star (playing 
the role of "secondary") are adjusted in a highly confined orbit (i.e., orbit with orbital separation 
less than a few solar radii). The secondary extends out to its critical Roche equipotential surface 
and small motions of thermal origin in its atmosphere lead to gas flow (through the inner Lagrangian 
point of the binary) into the potential well of the primary. This gas forms an accretion disk, which is 
an efficient machine at converting gravitational energy to radiation and heat, and at redistributing 
angular momentum. The accretion disk is disrupted at the Alfven radius and, eventually, the 
accretion flow is channelled towards the poles of the primary through magnetic columns (P94, § 
2.1). 

However, AE Aqr is characterized by some unusual properties. We mention here its violent 
flaring activity, with recurrence timescale ~ 1 h, from the ultraviolet (UV) (E&:94, §§ 2-4) to the 
radio frequencies ( Pastian, Dulk, &: Chanmugam 1988| , § 2). It is also remarkable that its spin 
period, Pnow = 33.08s (the shortest one among IPs), has been detected and verifled over a wide 
range of frequencies, from the optical ( Patterson 1979| , §§ 2, 6) to the 7-rays ( tVIeintjes et al. 1992| , 
§ 2, M&:94, §§ 2-3, 6). Further unusual properties of AE Aqr are discussed by, e.g., 195 (§ 1), 
WHG95 (§ 1), C&96 (§ 1), |Bruch fc Griitter 1997] , (§§ 2, 4-5), 197 (§ 1), K&97 (§ 1), and WKH97 
(§ 1). 

Concerning the subject of the present investigation, emphasis is given on the puzzling rapid 
spin-down time rate, Pnow = 5.64 x 10~^^ss~^, reported by dJ&:94 ( §§ 6-7, Table 4). The 
implications of such a spin-down have been discussed by several investigators (see, e.g., dJ&:94, § 
8; WHG95, § 1; [Beskrovnaya et al. 199^ , hereafter B&96, § 4; 197, §§ 1-3; 198, §§ 2-5). Especially 
in 197 (§ 1.2), the so-called "spin-down problem" is deflned clearly. In particular, when calculated 
classically, the observed spin-down time rate Pnow implies a spin-down power of the white dwarf 
(197, eq. [1]) 



l34, 



T ~ -6 X 10^''/5o-P-P"''ergs~' ~ -lO^^ergs 
where /50 is the moment of inertia of the white dwarf in units of 10^'^ 



(1) 

gcm^, and P and P are 
expressed in units of Pnow and Pnow, respectively. Comparing T with the observed luminosities 
(see, e.g., 197, Table 1), we realize that it exceeds the UV {Zuv) and X-ray (£x) luminosities of AE 
Aqr by a factor ^ 10^. This is in contradiction with the basic postulate of the IP model concerning 
the accretion nature of the primary's radiation (for details on this issue, see 195, § 2.3), and exactly 
this contradiction has been termed as "spin-down problem" . 



In the following two sections, we shall discuss some scenarios proposed for solving this problem; 
we shall also describe the scenario considered in the present investigation. 
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2. Dealing with the spin-down problem 

It is of interest, at this point, to mention some models concerning the magnetic, radiating, and 
accreting properties of AE Aqr; although some of these models are not directly connected to the 
spin-down problem, they help us to have a better understanding of the whole matter. In 195 (§§ 
2-3, Table 1), the state of the primary (under the basic assumption that it is a white dwarf) is 
taken to be simultaneously as ejector, propeller, and accretor. Then the spin-down power may be 
lost in the form of relativistic stellar wind ejected from the surface of the white dwarf (ejector) or 
transmitted to the accretion flow near the Alfven surface (combination of accretor and propeller; 
see also EH96, § 4.3). 

In K&97 (§ 8), it is concluded that if the rapidly rotating white dwarf in AE Aqr has a 
surface magnetic field of ~ 100 T, then magnetic pumping is a serious candidate for the particle 
acceleration needed to explain the observed radio outbursts. Magnetic pumping (K&:97, § 4 and 
references therein) is one of a number of stochastic particle acceleration mechanisms, known also as 
Fermi-type mechanisms, whose essential property is the combination of a particle interaction with 
a spatial or temporal change in the magnetic field — which causes the particle energy to change 
— with a randomizing process (see also EH96, § 5). 

WKH97 (§§ 2-3, 6) discuss in detail the so-called "diamagnetic blob model" (for some related 
topics of interest, see also [King fc Lasota 199l| , [Wynn &: King 19^ , [King 199"^ , |Wynn fc King 



1995 ). According to the authors, the discrepancy between the large spin-down power of the white 
dwarf and the observed luminosity of the system is explained by the kinetic energy carried away 
by the ejecta. This model does also predict a magnetic dissipation time rate ~ 10^^ ergs~^. There 
are some processes that may be triggered by such dissipation, but they could account for only 
^ 1 percent of the dissipated energy. A significant fraction of this energy may go into particle 
acceleration and production of 7-rays. 

In 197, a neutron star model is proposed for the primary in AE Aqr. The author concludes 
that, in the frame of the wind- fed accretor model, there are no particular problems with the 
interpretation of the observed spin-down time rate of the primary (remember that the moment of 
inertia of a neutron star is ~ 5 orders of magnitude less than that of a white dwarf). However, this 
scenario reveals a very intriguing problem (197, § 4): if the primary in AE Aqr is indeed a neutron 
star, why its appearance in the optical UV is very likely to that of an accreting white dwarf? 

In a subsequent investigation, Ikhsanov (198, §3) recognizes that, in spite of some progress, the 
origin of the observed spin-down of the primary in AE Aqr cannot be regarded as well understood so 
far; in this situation the investigation of alternative possibilities (if they exist) remains important. In 
this view, the author (198, §§ 4-5) discusses a pulsar-like white dwarf model in a binary system. In 
particular, he explores the possibility to solve the spin-down problem of AE Aqr within the "ejector 
approach" , assuming that the major part of the rotational kinetic energy of the white dwarf is spent 
in the generation of low frequency magnetodipole waves and particle acceleration. As the author 
discusses, the main consequence of this approach is a relatively strong magnetic field of the white 
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dwarf, which makes this model to differ radically from previous ones. The estimated magnetic 
moment (~ 10^^ G cm^) is ~ 2 orders of magnitude greater than those of previous models; and the 
question arising concerns the consistency of such an estimate with corresponding observations (see 
198, § 5 for interesting details and references therein; especially. Cropper 1986| , Stockman et al 



1992, B&96). 



Recently, Ikhsanov (199) describes a scenario according to which "magnetic amplification" can 
be achieved by the action of differential rotation, appearing, in turn, as consequence of gravitational 
waves (and resulting angular momentum loss); if so, then the magnetic field inside the star winds 
up, on a timescale of ~ 1 month, up to ~ 10^ -i- 10^ G and the star spins down predominantly due to 
magneto-dipole waves and particle acceleration. Furthermore, in the framework of the pulsar-like 



white dwarf model (as in 198, 199), Ikhsanov (2000; see §§ 3-4 for interesting details) suggests the 



existence of a dead disk around the white dwarf magnetosphere. The main feature that makes the 
suggested model different from those previously discussed, is the hot boundary layer at the inner 
disk radius. Due to this feature, the model sheds new light to some puzzling properties of AE Aqr. 

I ,3 

Meintjes & de Jager (2000), in order to explain the multiwavelength emission from AE Aqr, 

propose a model on the basis of the assumption that the compact component of this binary system 
consists of a white dwarf with a surface field length of ~ 1 MG, rotating through a clumpy ring 
near the circularization radius. After a detailed analysis, they conclude that the whole reservoir of 
the spin-down power has gone into the production of very high-energy particles and 7-ray emission. 



Finally, Choi &: Yi (2000) point out that the electromagnetic dipole radiation with surface 
magnetic field ^ 5 x 10^ G and the propeller action with accretion time rate ^ lO^^gs"^ are 
unable to account for the observed spin-down torque. So, the authors propose that the spin- 
down power does not transform into the observed electromagnetic radiation and they consider the 
gravitational radiation emission as an alternative spin-down mechanism. The authors remark that 
this mechanism could be an attractive one since the resulting spin-down power does not need to go 
into the observed electromagnetic radiation, which effectively avoids the long standing question of 
nondetection of the observed large spin-down power. 



3. The turn-over scenario 

In the present paper, we try to resolve the spin-down problem in terms of intrinsic prop- 
erties (like angular momentum, moments of inertia, differential rotation, and magnetic field) of 
the primary in AE Aqr. In particular, we propose an axisymmetric, differentially rotating, mag- 
netic, white-dwarf model with its surface magnetic field causing a slow "secular spin-down" (to 
be distinguished from the fast spin-down observed now) and, accordingly, a slow "secular angular 
momentum loss" . This model participates simultaneously in the accretional activities of the binary 
system AE Aqr; and, apparently, its secular activities are screened by its accretional activities 
as long as the former remain weak in comparison with the latter. Under certain circumstances. 
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however, such secular activities can cause certain quite intense phenomena. In this view, we shall 
try to document the aspect that the fast spin-down of the primary in AE Aqr, observed now, be- 
longs to such a category of phenomena. In particular, we shall give numerical results showing that 
such fast spin-down may be the result of a nonaxisymmetric transition from an initial axisymmet- 
ric differential rotation state to a terminal nonaxisymmetric rigid rotation state, due to turn-over 
of the magnetic symmetry axis with respect to the invariant angular momentum axis (the latter 
taken to be the z axis). Our results mainly show that a fast spin-down does not always imply a 
large spin-down power; instead, under certain circumstances as this of turning over, the spin-down 
power corresponding to a large spin-down time rate (as that observed in AE Aqr) is ~ 2 orders of 
magnitude less than that estimated classically. 



We shall use hereafter definitions and symbols identical to those in Geroyannis & Papasotiriou 



1999, hereafter GP99, Table 1) 



According to the so-called "turn-over scenario" (hereafter TOV scenario) , the model undergoes 
an axisymmetric "early evolutionary phase", during which the combined effect of rotation and 
(toroidal plus poloidal) magnetic field derives configurations with moment of inertia along the 
rotation axis /s (coinciding with the symmetry axis), 133, greater than the moment (s) of inertia 
along the other two principal axes 1\ and -^11 = -^22; i-e., ^33 > At some particular time, 
however, which signals the beginning of a "late evolutionary phase", secular angular momentum 
loss yields configurations with 133 (L^^) = l\\{Lx) (relation valid for a particular angular momentum 
Lx) and, after that particular time, configurations with /33(-L) < 1\\{V) when L < L^- Thus, a 
"dynamical asymmetry" is established in the sense that the moment of inertia along the rotation 
axis becomes less than the moment(s) of inertia along the other two principal axes. This inequality 
is due to the action of toroidal magnetic field, which tends to derive prolate configurations; while, 
on the other hand, rotation and poloidal magnetic field tend to derive oblate configurations. 

Eventually, secular decrease of rotation due to surface magnetic field gives rise to dynamical 
asymmetry, 1\\ > 133. However, dynamically asymmetric configurations tend to turn over sponta- 
neously, so that to subsequently rotate along axis with moment of inertia greater than 133. They 
become, therefore, "oblique rotators" ; the well-known "perpendicular rotator" is an oblique rotator 
with "turn-over angle" (or, equivalently, "obliquity angle": the angle of the magnetic symmetry 
axis with respect to the invariant angular momentum axis), X) equal to 90°. The reasons causing 
turn-over of dynamically asymmetric configurations have been studied by several investigators (see, 
especially, [Mestel fc Takhar 1972 , hereafter MT72); details on this issue are given in § 6. 



Within the framework of TOV scenario, we shall study the case in which the terminal nonax- 
isymmetric rigidly rotating model becomes perpendicular rotator on a "turn-over timescale", Itov^ 
comparable with the "period timescale", tp = Pnow/Pnow] i-e., tTov ~ tp- Although any inter- 
mediate oblique model with x < 90° undergoes a complicated rotation along axis passing through 
the center of mass but not coinciding with one of its principal axes (MT72, §§ 2-3), the terminal 
model rotates along its Ii principal axis, which axis is now coinciding with the invariant angular 
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momentum axis taken, in turn, to be the z axis. Apparently, its angular momentum vector remains 
equal to that of the starting model (MT72, §§ 1-2), while its angular velocity vector obtains same 
direction with that of the starting model (coinciding with the angular momentum direction) but 
different magnitude. 

During the turn-over phase, the turn-over angle increases spontaneously up to 90° on timescale 
t-roV) since the rotational kinetic energy of the model decreases from a higher level with turn-over 
angle x ^ 0° to a lower level with turn-over angle x — 90° • At this level, the model reaches the state 
of least energy consistent with its prescribed angular momentum and magnetic field. Accordingly, 
the excess energy due to the differential rotation field of motions, defined by the angular velocity 
component along the (spontaneously turning over) /a principal axis (coinciding, in turn, with 
the magnetic symmetry axis), is dissipated down to zero under the action of turbulent viscosity 
in the convective regions of the model (MT72, §§ 3-4). Thus, differential rotation along the I3 
principal axis decreases down to zero and all of its effects disappear. 

Furthermore, it seems difficult for an oblique model to sustain differential rotation along its Ii 
principal axis mainly due to the destructive action of the poloidal magnetic field (see, e.g., [Moss 



Mestel, fc Tayler 1990| , hereafter MMT90, §§ 1-3; also, [Moss 1992| , hereafter M92, §§ 1, 5-6). In 



simple words, in a model with magnetic field symmetry axis inclined with respect to the rotation 
axis, i.e., x > 0°, nonuniform rotation in the electrically conducting stellar material shears the 
field and generates magnetic torques. As a result, there is a competition between the efforts of 
the magnetic stresses to remove rotational nonuniformities, and those of the rotational velocities 
to bury and destroy magnetic flux. If, then, the magnetic field and the electrical conductivity 
have appropriate values (issues to be discussed in § 7; see, e.g., MMT90, § 4; also, M92, § 6), 
then the magnetic field prevails and removes all the nonuniformities of rotation. Consequently, 
the terminal model rotates rigidly along its Ii principal axis with angular velocity However, 
differential rotation behaves like a "rotational kinetic energy capacitor" and its decrease down to 
zero discharges all this excess energy (since the rigidly rotating perpendicular rotator has a much 
lower "rotational kinetic energy capacity"). Hence, differential rotation decrease leads to releasing 
the excess rotational energy on a "differential rotation decrease timescale", tz)_RD ~ ^TOVr Our 
task here is to verify that the corresponding energy release time rate is ~ 10^ times less than that 
estimated by equation (||). 

In detail, we assume that the model of TOV scenario undergoes a nonaxisymmetric transition 
from its axisymmetric DRS to its nonaxisymmetric perpendicular RRS. The starting model has 
mass MjjR = 0.89Mq, angular momentum Ljjji = L^x such that L(Pnow) < L^x < Lx, rotational 
kinetic energy Tdr = Txx, spin period Pdr = Pxx < Pnow, magnetic perturbation parameter 
and magnetic flux foR- The terminal model has same mass, Mrr = M^r, same angular 
momentum, Lrr = Lxx, and same magnetic flux, /rr = foR- If so, then the fitting between 
the starting model and the terminal model of such a "nonaxisymmetric DRS-to-RRS transition" 
(NXDRT) can proceed through the magnetic perturbation parameter Hrr of the terminal model. 
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The equations 

Lrr = Lxx = lRR,ll{hRR, Lxx)fli, (2) 

TRR{hRR) = -lRR,ii{hRR, Lxx)fll, (3) 

give, for various values Hrr of the magnetic perturbation parameter, the corresponding angular 
velocity component Jli along the Ii principal axis (with moment of inertia lRR,ii{hRR, L^x)] details 
on the calculation of this quantity are given in § 5) and the rotational kinetic energy Trr. If, 
furthermore, NXDRT takes place under constant magnetic perturbation parameter, Krr = h^R 
(which seems to be a plausible assumption), then the system of equations (§)-(^) gives the proper 
values for Trr, and the period, Prr = 27r/r2i, of the terminal model. Then an estimate of the 
timescale t_D_RD needed for the spin period to increase from the value Pxx to the value Prr is 

, Prr Pxx / .\ 

toRD ^ — 5 ; (4) 

and an estimate for the NXDRT power is 

rh '^xx — Trr DEdrd 

1 ~ = . (5j 

toRD toRD 

In view of TOV scenario, therefore, the spin-down problem is drastically simplified, since, as our 
numerical results reveal, the NXDRT power is ~ 10^ times less than the classically estimated 
spin-down power of equation (||). 

Finally, the following remark is of interest. In order to keep computations within the two- 
dimensional framework of the so-called "complex-plane iterative technique" (CIT) (GP99, §3 and 
references therein) , we simulate the terminal rigidly rotating models by a constant-mass sequence of 
magnetic models with zero rotation (i.e., with rotational perturbation parameter z; = 0) along the /s 
principal axis. Apparently, all the members of this nonrotating magnetic constant-mass sequence 
are prolate along the /s principal axis. Thus, our computations retain acceptable reliability, in 
describing the phenomenon under consideration, as far as the component $^3 of remains small 
with respect to the component f^i, the turn-over angle x remains near ^, and the angle 7 of the 
angular velocity direction with respect to the angular momentum direction remains small. For the 
perpendicular rotator we have ^3 = 0, x = 90°, and 7 = 0°; thus, the perpendicular rotator can 
be described reliably by the simulation adopted. 



Details on the overall numerical method 



The axisymmetric model, initially used in TOV scenario, is computed by (the two-dimensional) 
CIT. In a previous investigation ( |Geroyannis"T99T| , hereafter G91, §3), CIT has been developed and 



used for computing rotating visco-polytropic models; this term indicates that the differential rota- 
tion of such models results from the particular equations governing the viscosity of their material. 
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CIT has been extended for computing rotating white dwarfs ( [Geroyannis fc Papasotiriou 1994 , here- 



after GP94), and rotating white dwarfs with toroidal magnetic field ( peroyannis &: Papasotiriou 



1997| , hereafter GP97). Recently (GP99), CIT has been further extended for computing rotating 



white dwarfs distorted by both toroidal and poloidal magnetic fields. An almost complete set of 
quantities involved in the computations is given in GP99 (Table 1). As we have already clarified, 
in the present paper we shall use definitions and symbols identical to those in GP99. 

The algorithm of CIT is analyzed in G91 (§ 3; see also GP99, Fig. 1). The structure and 
the rotational behavior of the axisymmetric model is identical to that of GP94 (§2). The toroidal 
magnetic field, Ht = Htu^, is defined as in GP97 (§2, eqs. [2.1]-[2.2]; see, also, GP99, eq. [2.7]); 
the poloidal magnetic field, Hp = HpT^u.-^ + HpzUz, is defined as in GP99 (§2, eqs. [2.8]-[2.9]; w, 
(f, and z are the cylindrical coordinates). On the base of these definitions and, also, of the units 
and the corresponding dimensionless quantities defined in GP99 (Table 1), some straightforward 
algebra shows that 

3 

Ht = Htu^ = {kwQ)u^ = [2TrVGaBy^]Vhs (^^ - \) ' u^, (6) 

/ k dg\ / k k dg\ ,_, 

Hp = Hp^u^ + HpzUz = (^-^-^j + (2-, + -w—j u., (7) 



Hp^ = [27rVGaByl]^^^ l-^s^J^^ - , (8) 



Hp. = [2.VGaByl]J^^ U - ^) % 3.<1> ) , (9) 



yj V yo 

where s = w/a and Q = z/a are the dimensionless cylindrical coordinates (a is the well-known 
white-dwarf length unit). For configurations declining slightly from sphericity, we can use the 
approximative expressions 

^ ^ VT^^, (10) 

OS 

where ^ = r ja and v = cos are the dimensionless modified spherical coordinates (taken together 
with the angle (p; on the other hand, r, and (p are the usual spherical coordinates). 

The average surface poloidal magnetic field 



^ {HpiCs,!^)), = {^Hl^{is,y) + Hl,{is,u))^ , (12) 

where the subscript v denotes the average over this coordinate, is somewhat prescribed to the value 
B'P by either the theory or the observations. So, there is need to define a "surface zone" with "base" 
at (dimensionless) radius ^s, say. For white-dwarf models, a common practice is to consider that 
the points below the boundary where electrons become degenerate — i.e., the "transition layer" 
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— are located at distance from the center (see, e.g., ^hapiro Sz Teukolsky 1983 , hereafter ST83, 
§4.1). So, the transition layer coincides with the base of the surface zone. The temperature at 
this transition layer (which is in turn equal to the interior temperature), 1^, is given in terms of 
the luminosity, £, the mass, M/Mq, and the composition, X, Y, Z, appearing either explicitly or 
through the quantities 

m = 5 1 — , (13) 

= x + iy + iz' 

where m is the "mean molecular weight" and rrie the "mean molecular weight per free electron", 
both measured in proton masses (see, e.g., ST83, §2.3). In particular, is given by the relation 
(ST83, eq. [4.1.11]) 



Z{l+X)mi£ 



1 

3.5 



- I 5.7 X lO^m J| ' ' ^^^^ 
with corresponding density at the transition layer (ST83, eq. [4.1.10]) 

Qs = 2.4 X 10~VeTj-^ (16) 
Thus, the radius results as root of the nonlinear algebraic equation 

QiCs) -Qs = 0, (17) 

or, equivalently. 



1 \ 

^ (6) - - -Qs = 0, (18) 
which can be solved by a standard algorithm of numerical analysis. 



The fitting between the prescribed average surface poloidal magnetic field, B^, and the "proper 
poloidal magnetic parameter", P^, can proceed through the nonlinear algebraic equation 



= {^H^^ih; 6, i^) + Hj,{h- I3l- 6, v)) ^ , (19) 

which, for prescribed /i, can be solved by a standard algorithm of numerical analysis and give /3j 
as root. 

Concerning the computation of constant-mass sequences, the white dwarf case is quite different 
to the polytropic case. In the polytropic case, one can construct a constant-mass sequence simply by 
keeping a particular fixed value for the fundamental model parameter n, i.e., for the polytropic index 
(see, e.g., ^hapiro fc Teukolsky 1990 , §2; also, Galli 1992, §3). On the other hand, a constant-mass 
sequence of white dwarfs consists of models with properly varying fundamental model parameter 
2/0 or, equivalently, This is due to the fact that the central density Qc depends explicitly on yo 
through the well-known relation 

1 \3/2 



Byt ( 1 - ^ ) ■ (20) 
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Thus, a rotating magnetic white dwarf may attain a specific mass only for a specific combination 
of the fundamental model parameter yo and the rotational perturbation parameter v. 

It is worth emphasizing that white dwarfs with given 4j- may not be able to attain a specific 
mass at all, since there are lower and upper mass limits; namely, the lower mass is that of a 
nonrotating nonmagnetic model, and the upper mass (for given magnetic perturbation parameter 
h) is the mass of a critically rotating model. 

In order to compute constant-mass sequences of white dwarfs, we use CIT in an iterative 
algorithm, called "search for specific mass" (SSM). This algorithm finds the value VM,h for which 
a model with given and h has a specific mass M/Mq. Alternatively, SSM can find the value 

hM,v=o for which a nonrotating model (i.e., i; = 0) with given 4j has a specific mass M/Mq. Using 
SSM we can compute several constant-mass rotating sequences of white dwarfs. By calculating the 
angular velocity and the angular momentum of each member of a rotating sequence, we can plot 
the curve i}{L) corresponding to a specific mass. However, in order to construct "complete curves", 
we need to know the extrema of the fundamental model parameter ^ , for which the sequence can 
still attain the specific mass. The maximum value {^)max corresponds to the rightmost point of 
the curve ^{L), that is, to an almost critically rotating model. A model with slightly larger ^ will 
fail to attain the desired mass, since its mass is less than the desired one even in critical rotation. 
On the other hand, the minimum value (^)min corresponds to the leftmost point of the curve, 
that is, to a model with rotational kinetic energy far from its critical value. A model with slightly 
smaller 4^ cannot attain the desired mass, since its lowest mass is greater than the desired one. 

To calculate the extrema of ^ we have developed the so-called "search for extrema of the 
basic model parameter" (SEB). Starting from an intial guess for each extremum, this algorithm 
uses SSM iteratively to calculate the corresponding extrema. Alternatively, SEB can find the 
corresponding extrema of ^ for a constant-mass nonrotating sequence, by properly adjusting the 
magnetic perturbation parameter h. 

To summarize, the computation of a rotating magnetic constant-mass sequence of white dwarfs 
consists of the following procedure. First, we specify certain values for the mass M/Mq, the 
magnetic perturbation parameter h, and the differential rotation reduction factor Fr- Alternatively, 
we take v = — nonrotating magnetic constant-mass sequence — and we use in the fitting 
procedure a varying h. Second, we use SEB to calculate the corresponding (\)min and (\)max- 
Third, we use SSM to compute several models within the already specified range of values. By 
these results we can study the behavior of the corresponding curve fl{L). Alternatively, we can 
study the full range of values {^,h) for which a nonrotating magnetic sequence retains constant 
mass. 
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5. Calculation of the moment of inertia of the perpendicular rotator along its 

rotation axis 

The quantity lRR,ii{hRR,Lxx)-, involved in equations denotes the moment of inertia 

of the perpendicular rotator along its rotation axis Ii. The problem of calculating this moment is 
three-dimensional, since the equilibrium structure of the perpendicular rotator is represented by 
an ellipsoidal figure. Thus, a systematic numerical treatment of this subject consists in developing 
an iterative three-dimensional method, like that of Hachisu (1986)| , for computing equilibrium 



structures distorted by both rotation and magnetic field. Our plan (for the near future) is to 
modify and extend CIT in order to obtain the "three-dimensional CIT" (3DCIT), able to make the 
job required. At present, however, 3DCIT is not available. So, there is need to obtain a reliable 
estimate oi lRR^ii{hRR, L^x) by implementation of the existing "two-dimensional CIT" (2DCIT). 

We first clarify that 2DCIT can calculate (within a prescribed numerical accuracy) the mo- 
ments Iii{hRR) and /33(/i_Ri?) of a nonrotating magnetic constant-mass sequence, and the moment 
lRR,33ihRR, Lxx) of an axisymmetric rigidly rotating magnetic constant-mass sequence. Indeed, in 
both cases the corresponding equilibrium structures are spheroidal figures; thus, they are involved 
in a two-dimensional problem easily handled by 2DCIT. 

Now, there is still a chance to obtain a reliable estimate of the so-called "rigid rotation ampli- 
fication ratio", Ar{Ii), resulting from the definition 

^r\i-i) — P77 X l^-Lj 

^uyriRR) 

(with the substitutions i = 1, n = 11), within the framework of 2DCIT. It is expected that 
^T-(-^i) > 1) since rotation does prolate further the configuration of the perpendicular rotator along 
the /a principal axis and, so, does increase the moment of inertia along the rotation axis Ii from 
the value /ii(/ir_r) (magnetic field only) to the value lRR,ii{hRR, Lxx) (magnetic field plus rigid 
rotation). To start with, let us consider our problem in view of a perturbation method, like that 



of Geroyannis &: Hadjopoulos (1990, hereafter GH90) . In this view, the moments lRR,ii{hRR, L, 



and lRR,33{hRR, Lxx) are given by the expansions (see GII90, eqs. [3.6]-[3.8], [4.1], and other similar 
expansions therein) 

lRR,ii{hRR, Lxx) = liiisph) + AIii{mag) + AIii{rot) + AIii{mag * rot), (22) 

where the argument sph shows the spherical (i.e., the nondistorted) contribution, the argument 
mag gives distortion due to magnetic field, the argument rot represents distortion due to rotation, 
and the argument mag * rot shows mixed distorting action of both magnetic field and rotation. 
The significant point is that, when approximating the problem by a first-order perturbation method 



(see, e.g., Geroyannis & Hadjopoulos 1992, §§ 1-2), then, apart from the spherical contribution, 
we take into account distortions due to magnetic field of first order in the magnetic perturbation 
parameter h, 

AIii{mag) ~ hM,i, (23) 
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and distortions due to rotation of first order in the rotational perturbation parameter v, 

M,i{rot) ~ vRu, (24) 

while distortions of order equal to or higher than /i^, u^, hv (mixed distorting terms), are omitted. 
It appears therefore that, within the framework of a first-order approximation, the distortions due 
to magnetic field and rotation become decoupled and the corresponding calculations 

lRR,ii{hRR, L^x) ^ hi{sph) + Mii{mag) + Mii{rot) (25) 

can be undertaken in separate steps concerning sph-terms, mag-terms, and rot-terms. Furthermore, 
combination of equation (21), defining Ar{Ii), and of the approximating equation ( [25| ) gives 

^ ^ liijsph) + AIii{mag) + AIii{rot) Alajrot) 
Iii{sph) + AIii{mag) hiihRR) ' 

which says that either Ar{Ii) is known and then AIii{rot) can be calculated by 

AUrot) ^ -.\ (27) 



or AIii{rot) is known and then Ar{Ii) can be calculated by equation (^61). 

Now, there is an obvious "invariance" concerning rotational distortion with respect to the 
principal axes Ii, I2, I3, when departing from an initial spherical (i.e., nondistorted) configuration. 
In particular (and for same rotational perturbation parameter v), when taking I3 as rotation axis, 
we calculate a distortion A/33 (rot) equal to the distortion Al22{rot) calculated when taking I2 as 
rotation axis, and, in turn, equal to the distortion AIii(rot) calculated when taking Ii as rotation 
axis. Because of this invariance, the distortion AIii{rot) is equal to the distortion A/33 (rot) 
calculated by equation ( p7| ) (with the substitutions i = 3 and ii = 33). We remark that all the 
quantities involved in equations (|26|)-(p7|) can be calculated within the framework of 2DCIT; so, 
the moment of inertia /r_r,ii(/i_r_r, /^xx) = ^r(/i)/ii(/i_R_R) can be eventually calculated by 2DCIT. 

The next significant point is that the angular velocity 0,i of the perpendicular rotator, 

= A.{h)Iu{hRRy ^^^^ 

should give rotational perturbation parameter v, 

- 2nGByl{hnR.) ^^^^ 

(where the function i/q = yo{h) is known from our nonrotating magnetic constant-mass sequence), 
equal to the value of v used in equation (p7|) (with the substitutions i = 3 and ii = 33) for evaluating 
A/33 (rot), which is taken, in turn, from the known function v = v{L) holding for our axisymmetric 
rigidly rotating magnetic constant-mass sequence. Accordingly, a fitting is required between the 
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value of V resulting from the approximation (25)-(2£) and the value v = v{L) used for deriving this 



approximation. A simplified iterative algorithm undertaking this fitting can be constructed on the 
basis of the remark that the so-called "rotational distortion coefficient" Ra in equation ( |2^ ) varies 
slow with the angular momentum, while v varies fast. In other words, Ra can be assumed constant 
over a narrow range [L{Pnow),Lx\ of angular momenta. In detail, this "simplified algorithm" (S) 
has as follows. 

Step 1(S). Take L = Lrh = L^x for the angular momentum and h = Hrr = h£)R for the 
magnetic perturbation parameter; then calculate the rotational perturbation parameter v^^^ from 
the known function v = v{L) holding for our axisymmetric rigidly rotating magnetic constant-mass 
sequence, v^^^ = v{Lxx)- 

Step 2(S). Calculate the rigid rotation amplification ratio Ar{Iz) from the relation (eq. |^l[ 
with i = 3 and ii = 33) 

. lRR,33{hRR,Lxx) .„„x 
Arih) = J— 77 N , l3Uj 



the distortion A/33 (rot) from the relation (eq. with i = 3 and ii = 33) 



A/33(™t) \ (31) 



and the rotational distortion coefficients Ru, Rs3 from the relation (eq. [gj] with ii = 33) 

-Rii = R33 = — ^^Y) — ^' ^^"^^ 

Step 3(S). Calculate the distortion A/ii(rot)^^) from the relation (eq. (2j] with ii = 11) 

A/11 (rot) W = iiiiuW, (33) 
the rigid rotation amplification ratio Ar{Ii)^^^ from the relation (eq. with i = 1 and ii = 11) 

+ (34) 

and take lRR^ii{hRR, Lxx)^^^ = Ar{Ii)^^^ Iii{hjiR); then calculate the angular velocity Q^l^ from the 
relation (eq. [p8[) 

^S'^ = 7 1(1) ' (3^) 

^RR,ll{ilRR, J-'xx)^ ' 
and the rotational perturbation parameter v^'^'^ from the relation (eq. |]29|| ) 

(2) _ {^h ) /ogx 

27,GByl{hRny ^^^^ 
where the function y = yo{h) is known from our nonrotating magnetic constant-mass sequence. 
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Step 4(S). Check if \ v^'^^ —v^^^ \ < tol, where tol is a proper tolerance. If yes, then iRR^uihjm, L^x)^^"^ 
is a satisfactory approximation of -^/j_r,,ii(/ihj?) L^x)- If no, then use the value v = v^"^^ for the next 
iteration, starting from Step 3(S). Remark that the iteration counter obtains now the value (2) for 
Step 3(S), where, however, the new estimate for v takes as iteration counter the value (3), So, 
when arriving again at Step 4(S), the comparison is made between the values v^^^ and v^^^; etc. 

The "complete algorithm" (C), on the other hand, takes into account even the slow variation 
of the rotational distortion coefficient Ra. This algorithm has as follows. 

Step 1(C). Take L^^-* = Lrr = Lxx for the angular momentum and h = Hrr = h^Ft for the 
magnetic perturbation parameter; then calculate the rotational perturbation parameter v^^^ from 
the known function v = v{L) holding for our axisymmetric rigidly rotating magnetic constant-mass 
sequence, v^^^ = v{L^^^). 

Step 2(C). Calculate the rigid rotation amplification ratio Ar{l3)^^^ from the relation 

^ri^sr = — riT. — ^ — ' 

33 ("-iJijj 

and the distortion Al3s{rot)^^^ from equation (p7|); then take Alii(rot)^^^ = A/33(rof)^^^ and 
calculate the rigid rotation amplification ratio Ar{Ii)^^^ from the relation 

+ (38) 

Step 3(C). Take lRR,n{hRR, L^^^) = Ar{IiY^^ In{hRR); then calculate the angular velocity Q,^^^ 
from the relation 

' ~ lRR,u{hRR,Lmy ^''^ 
and the rotational perturbation parameter v^"^^ from the relation 

y(^) = ^li^J (40) 

27rGByUhRRy ^^^^ 
where the function i/q = yo{h) is known from our nonrotating magnetic constant-mass sequence. 

Step 4(C). Check if |f (2)— f (^)| < tol, where tol is a proper tolerance. If yes, then lRR^ii{hRR, L^^^) 
is a satisfactory approximation of lRR,ii{hRR, Lxx)- If no, then calculate the root L^^^ of the alge- 
braic equation v{L) = v^'^\ and use the value L = L^^^ for the next iteration, starting from Step 
2(C). Remark that the iteration counter obtains now the value (2) for Steps 2(C) and 3(C), where, 
however, the new estimate for v takes as iteration counter the value (3), v^'^\ So, when arriving 
again at Step 4(C), the comparison is made between the values f^^-* and t;^^); etc. 



6. Turn-over and differential rotation decrease in nonaxisymmetric transitions 



MHD studies of rotating stars have been undertaken by several investigators (see, e.g., Chan 



drasekhar 1956a; Chandrasekhar 1956b; Chandrasekhar & Kendall 1957; Weber & Davis 1967; 


Parker 1970 


; Chanmugam & Gabriel 1972; MT72; 


Monaghan 1973; Wright 1973; Wright 1974; 


Mestel & Moss 1977; Moss 1977a, hereafter M77; 


Moss 1977b; Chanmugam 1979; Dicke 1979; 


Moss 1979; 


Moss 1981; Galea & Wood 1985; Pitts k Tayler 1985; Talmadge, Richter, & Fischbach 



1985; Tassoul Sz Tassoul 1986a, hereafter TT86; Tassoul Sz Tassoul 1986b; Mestel fc Weiss 1987 



Mestel, Moss, fc Tayler 1988| , hereafter MMT88; [Tassoul fc Tassoul 1989| , hereafter TT89; MMT90; 
Charbonneau fc MacGregor 199% hereafter CM92; M92; pharbonneau fc MacGregor 1993| ; [Mestel 
fc Shibata 199^ ; [Moss et al. T995| ; [Tavakol et al. 1995| ). 



In case of a model undergoing "axisymmetric DRS-to-RRS transition" (XDRT), the effect of 
mutual interaction of mixed magnetic field and rotation has become an issue of long debate among 
investigators asserting different aspects. In particular, a central question under debate has been 
the following: is there, in general, an unobtrusive magnetic field that can enforce uniform rotation 
and laminar meridional streaming? First, the study by TT86 (§ 8) gives negative answer to this 
question, unless turbulent viscosity is also involved in the process. Second, however, MMT88 (§ 
1, 4) argue that differential rotation — unless it is such that all parts of any poloidal field line 
have same angular velocity (situation known as "isorotation" ) — will produce a toroidal field from 
the poloidal field. So, if differential rotation is maintained, the toroidal field will build up steadily 
until either the field is able to diffuse through the matter or the Lorentz force becomes so strong 
that it reacts back on differential rotation. Eventually, unless the magnetic field is very weak, the 
influence of the Lorentz force will be so strong that only very slight differential rotation will be 
possible. Third, in turn, TT89 (§ 4) reconsider this matter and find that, although a fossil magnetic 
field may sometimes enforce almost uniform rotation at some places, it never acts uniformly in the 
whole stellar radiative envelope; instead, it always tends to generate large nonuniformities in both 
the angular velocity and the toroidal magnetic field — unless a large eddy viscosity is assumed to 
participate in the process, which then counters the effects of the magnetic field. The apparently 
conflicting results of TT86 and TT89, on the one hand, and MMT88 (see also MMT90, § 4), 
on the other, have been discussed in detail by CM92 (§ 4.4). The authors remark that it is the 
anchoring of all poloidal field lines in a rigidly rotating core that is ultimately responsible for 
the establishment of a state of near solid body rotation throughout the envelope. Therefore, the 
extend to which a weak poloidal field can maintain or enforce solid body rotation throughout the 
envelope depends critically on assumptions made regarding the behavior of the poloidal field at the 
boundaries (i.e., anchored field vs. unanchored field). However, they do also remark that even for 
fully core-anchored poloidal fields, quasi-steady differentially rotating configurations can sometimes 
materialize and survive for significant time intervals. 



In the present investigation, we focus our attention on NXDRTs only, mainly for three reasons; 
namely, when considering XDRTs (1) it is not clear how much strong should the magnetic field be 
so as to remove diff'erential rotation, (2) it is not clear when and why do such transitions start, and 
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(3) there should be strict axisymmetry in such models, since otherwise they could evolve to oblique 
rotators. 

To begin with, an approximative estimate for the turn-over timescale, tTOV, is given by (MT72, 
eq. [54] ; see also some interesting remarks regarding differential rotation at the last two paragraphs 
of § 3) 

where the subscript sz denotes that the corresponding average is calculated over the surface zone 
with base at the transition layer, and top at the boundary, ^i, of the star; averages without 
subscripts are calculated over the whole star, i.e., they are global averages. This estimate results 
as the mean energy per unit volume available for dissipation, DEuRo/y {DEdru is the excess 
rotational kinetic energy due to differential rotation, as in eq. |p, and V is the volume of the 
configuration), divided by the mean energy dissipation during the turn-over due to turbulent vis- 
cosity in the convective surface zone per unit volume and per unit time (i.e., the product of the 



last 5 terms in the right-hand side of eq. |41]; details on the derivation of such relation are given 
in MT72, eqs. [7], [49]-[50]). For the angular velocity, r23(s(^, z^)), we adopt the improved steady 
state Clement's differential rotation model (angular velocity constant on cylinders; as in GP99, eq. 
[2.3]) 



^3(5) = 03cw(s) = r^. 



3c 1 



. Y,a,e-Prb.s- , (42) 
\ 1=1 

where VL^c is the central angular velocity along the principal axis, u;(s) is the differential rotation 
function, ai and hi are the so-called "nonuniformity parameters" (their values for some typical cases 



are given in Geroyannis &: Hadjopoulos 1991 , Table 1), and F^. is the so-called "reduction factor" 
(note that lo = 1 when F^ = Q [rigid rotation]; and u: = minimum when Fy = 1 [situation known 
as "physical differential rotation"]). For the magnetic field, H(^,i^) = Hj(^,z^) -|- Hp(^,i/), we use 
equations (^)~(|9|). The remaining quantities in equation (^) are calculated under the assumption 
of slightly declining from sphericity (assumption valid when calculating Itov only). Accordingly, 
the density is approximated by q{^, u) ~ q{Oi the pressure is approximated by P(^, u) ~ P{i), the 



local sound speed, a, is approximated by (see, e.g., Musielak fc Fontenla 1989, § 2. a 



a{i,v) a{i) = , (43) 

the turbulent mixing length, A^, results from the "mixing-length theory" (MLT) — see, e.g., |Canuto| 



& Mazzitelli 1991, hereafter CM91, for interesting issues concerning MLT — as A^ = Aphp (CM91, 
§§ 2.8, 3.2), where the fine-tuning parameter Ap is taken to be of order unity, Ap ~ 1, and hp is 
the pressure scaleheight, 

A,«..)^A.«)^M«^^^|«); (44) 
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M(^) is the mass inside a sphere of radius ^, and the turbulent velocity, Vt, is approximated by 
(see, e.g., Brickhill 1990, hereafter B90, eq. [5]; in combination with CM91, eq. [38]) 



^*(e,^)=^^*(0=^^^- (45) 

It is worth remarking that, when adopting such approximations for hp and Vt, their average values 
over the surface zone for a white dwarf with mass M ~ 0.9Mq, radius R = a^i ~ 7 x 10*^ cm, and 
temperature ~ 4 x 10'' K are, respectively, 

{hp)s. - 3 X 10^ cm ~ ^ , (46) 

{yt)sz ^ % A ^ 5 X lO^cms-i (47) 

(compare with corresponding values given by |Fujimoto 1988| , § 2.1; see also [Fujimoto 1993| , § 4.1; 
and B90, § 4). 

Our numerical results show that, for model parameters chosen as above, the turn-over timescale 
is of order tTov ~ lO"^ Y- Thus, the magnetic poles move towards the equator (on trajectories 
coinciding with circumference quadrants, 27ri?/4) with an average speed, Vmp, given by 

Vmp = ~ IQ-^cms-^ (48) 

trov 

A question arising here is what kind of "globally acting" (i.e., acting over the whole star) "turn- 
over viscosity", utov, should substitute the "locally acting" (i.e., acting over the star's surface zone) 
"turbulent viscosity", fit, so as to induce energy dissipated over the time Itov-, DEtov, equal to 
the difference DEofid = T^x — Trr of the rotational kinetic energies of the starting and the terminal 
models. The reason for seeking a global viscosity, instead of the local turbulent viscosity, is that 
differential rotation, which causes this kind of viscosity, is a global property of the model. 

A first idea is to substitute the local turbulent velocity field, Vi, by a global "turn-over velocity 
field", Vtov ~ R/i-TOV-, which is of same order of magnitude with Vmp-, Vrov ~ Vmp- Then the 
turbulent viscosity (see, e.g., B90, eq. [8]) 

^t(e) = iQimio^phpic) ^ ^gio iv)^, {hp)^^ (49) 

(with = 1 as above) can be substituted by the turn-over viscosity 

f^ToviO ^ ^siO (Vtov) A^^ov (hp) ^ ^eiOj^^, (50) 

(since for our model {hp) ~ i?/10^) where Aj^py is a fine-tuning parameter. Furthermore, the 
energy dissipated due to such turn-over viscous friction per unit time, Dtov-, is given by (see, e.g.. 
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Geroyannis fc Sidiras 1995 , eq. [7]) 



Dtov = ^T^o?^lc I I ' /^Toy(0 - f ^) d^^^^^, (51) 

and assuming that Dtov remains constant over the time Itov the energy, DEtov^ dissipated due 
to such viscous friction over the time ixov is equal to tTovDrov ■ Typical values for our model 
are DEdjid ~ 10^'^ erg and, when taking Aj^Qy = 1, DEtov ~ 10^^ erg; thus, the required equality 
DEtov = DEdrd gives Al^Qy ~ 10 as proper order of magnitude for this fine-tuning parameter. 

A second idea is to develop a phenomenological model involving further details on the time 
evolution of the turn-over process (for similar considerations see M77, mainly §§3-4; see also p osa 



1990 , especially §§ 5-6, and Figs. 9-10). In the so-called "mixer-mixture model" (MMM), we 
assume that the turning over configuration acts as a "mixer" on the "mixture" configuration, 
which remains almost axisymmetric with respect to the invariant angular momentum axis, i.e., 
the z axis (the meaning of "almost axisymmetric configuration" is that the angle 7 between its 
angular velocity and its invariant angular momentum remains small: 7 < 1°, say; apparently, 7 
becomes zero after the termination of the turn-over; in the following, we shall assume that such 
a deviation of the mixture configuration from axisymmetry is negligible and we shall treat it as 
being axisymmetric with respect to the z axis). The mixer consumes its excess rotational kinetic 
energy due to differential rotation in the process of homogenizing the rotation of the mixture. In 
particular, when viewed in the mixture's system of reference, the mixer generates a velocity field, 
which carries mass elements across cylinders of different angular velocities and momenta. This is 
due to the fact that the rotational circular orbit of a mass element of the mixer, as seen in the 
mixer's system of reference, crosses several cylinders with different angular velocities, as seen in the 
mixture's system of reference. Thus, there is an interchange of angular momentum, which leads to 
uniform rotation of the mixture. 

For (1) the turn-over angle x (with = cos(x))j (2) the central angular velocity of 
the mixer, and (3) the differential rotation reduction factor F^- of the mixer, we adopt the time 
dependent scheme 

/ i + e(''-i)^ \ /l + e(''-i)^\ 

QiQinut) = L 1 QoiQini) - ( i ^ ^^5 I Qo(Qim), (52) 

where the so-called "threshold parameter", k, is taken equal to k = cos(45°) ~ 0.707, the so-called 
"delta parameter", 8, is chosen equal to 5 ~ 5, and the auxiliary function viit) varies linearly with 
the time in the interval [0, ^TOy]) 

if t < txov 



if t > txov 



When 5 is relatively large (5 = 50, say), the second term in the right-hand side of equation ( p2[ ) 
becomes negligible with respect to its counterpart, and QoiQini) can be simply set equal to the 
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initial value(s) Qini, i^xiim] ^ 1' ^3c[ini] = ^{L^^), and F^iini] = 1.00, or 0.50. However, as our 
numerical results show, the present study requires 5 ~ 5; since then the second term in the right- 
hand side of equation (52) cannot be neglected, QoiQini) must be normalized to the value 

Qini 
7 T+e("-i)'^ \ 

SO that Q{Qini,0) = Qini- 



QoiQ^m) = , (54) 



MT72 (§§ 2-3) use a linear perturbation analysis to describe the field of motions, ^ (MT72, 
eq. [44]), due to turn-over. In the framework of MMM, instead, we use a simplifying description 
according to which the turn-over speed at a particular point is proportional to the relative difference 
of the linear speeds due to differential rotation, Vmxr{£,j'^mxr,t) and Vmxti(,,i^mxt,t), in the mixer 
and in the mixture (where the rotation of the latter results as superposition of a progressively 
established rigid rotation about the z axis with angular velocity i^z [which is not involved in the 
calculation of the turn-over speed] and of a differential rotation about the z axis identical to the 
differential rotation of the mixer about the Is principal axis; for simplicity, we shall use the same 
symbols for the central angular velocity risc and the angular velocity field 1^3 for the differential 
rotations of both the mixer and the mixture, keeping in mind that the only difference concerns 
their rotation axes). 



Vtov = Vr 



mp 



Vmp 



Vmxt t^mxt, t) 
W{s{i, Vmxt), Fr{t))s{i, Umxt) - ^{sjj, l^mxr), Fr{t))s{(,, Vn 



(55) 



where Vmp is the expected order of magnitude for Vtov- Regarding the relation between the 
coordinates I'mxt and I'mxr of a particular point, some elementary algebra (where the axisymmetry 
of the mixture configuration with respect to the z axis is taken into account) shows that 



/ g \ _ I X + ^mxr if X + '&mxr < 90° 

ifmxt[i^mxr,X) - iggo _ + if ^ + > QQO ' ^^^^ 



and, accordingly, 




/ ^, ) V ^ " ^I^Xf — 90 l't^'7\ 

Assuming for simplicity that the coordinate v (without subscript) refers to the mixer configuration, 
we can write for the turn-over viscosity 

f^Tovit t) ~ ^QiOA^TOvyTOviC, t)A\ovhp{i) ■ (58) 
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Then the energy dissipated over the time tTOV, DEtov , is given by 

rtrov 



rtTOV 

DEtov = / DTov{t) dt = 47ra'^ x 

JO 







trov 



\J0 Jo V OS J s=s{i,u) J 



(59) 



remark that DTov{t) is now a function of time, while in the foregoing interpretation it was taken 
approximately constant over tTOV- Setting Al^Qy ~ 10 as above, and AJ^Qy = 1, we find DEdrt, ~ 
10^"^ erg and DEtov ~ 5 x 10^^ erg; thus, the required equality DEtov = DEdrd yields ^rov ~ 
2 X 10 as proper order of magnitude for this fine-tuning parameter. 

Eventually, we can adopt the above estimated values for the fine-tuning parameters, Ap = 1, 
AipQY = 10, Aj,Qy = 2 X 10, and try to find "optimal values" for the model parameters Lxx (i-e., 
the invariant value of angular momentum under which turn-over takes place), Bg (i.e., the average 
surface poloidal magnetic field), and 5 (i.e., the delta parameter of the time dependent scheme 
32|]). These optimal values result as the roots of the nonlinear system of equations 



^DRD = tTOV i-e-, equality imposed on relations (4) and 

DEdjit) = DEtov i-e-, equality imposed on relations (^) and (|59|), (60) 
DEDRD/toRD = DTov{t now ) ) 

which can be solved by a standard algorithm of numerical analysis. The last equation of this system 
has the physical meaning that W6 observe at the present time t^iow 

(such that P{tnow) = Pnow) an 

energy dissipation time rate, DTovitnow), more or less equal to its average value given by equation 
dll); in other words, we assume that the turn-over observed now is smooth and described by its 
average values. 

Finally, it is worth mentioning that the energy dissipation time rate due to the electron and 
ion viscosities of degenerate matter (see, e.g., Durisen 1973| , § 3, eqs. [4], [6], respectively), Deiv 



(calculated by eq. |^ with fiTOV substituted by the electron-ion viscosity ^e+i), is ~ 2 orders of 
magnitude less than Dtov] so, the former is totally screened by the latter and becomes insignificant 
in the framework of TOV scenario. 



7. Is the magnetic field strong enough to remove differential rotation? 

A condition stating when the magnetic field dominates over the differential rotation, and 
eventually removes it, has been set up for a perpendicular rotator by MMT90 (§§ 3-4; see also 
M92, §6). This condition involves the magnetic Reynolds number pi (MMT90, eq. [3.8a]), 
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and a measure p2 of the ratio V^^/V^qj^ (where Vas = ^sl \/'^^^s is the well-known Alfven speed 
at the base of the surface zone [see, e.g., CM92, eq. [13a]] and Vrot is the linear velocity due to 
rotation). 



P2 



(62) 



In these definitions, Urn is the coefficient of magnetic diffusivity (see, e.g., Tassoul 1978 , hereafter 
T78, Chapt. 3, eq. [125]; see also ST83, eq. [7.1.8]), 



47r(T ' 



(63) 



where a is the coefficient of electrical conductivity (see, e.g., T78, Chapt. 3, eq. [120]; see also 
ST83, eq. [7.1.6]); is taken here as the average value 17 = {Vtc + f7e)/2 = ilc(l + We)/2, where 
iVe = u){s = ^e) is the differential rotation at the equatorial surface; Bg is the average surface 
poloidal magnetic field (as in eq. []l^]); Qs is the density at the transition layer (as in eq. [p^); 
and R is taken here as the average value R = {Re + Rp)/^ = CK(Ce + Cp)/2, where ^e, ^-^e the 
dimensionless equatorial and polar radii, respectively. 



Using an argument set up earlier by Radler (1986) , MMT90 show that, in order for the magnetic 
field to dominate over the differential rotation in case of a perpendicular rotator, p2 must satisfy 
the condition (MMT90, eq. [4.1]) 



P2 > P2c ^ 



2/3 
Pi 



4/3 



(64) 



where p2c denotes a critical value for p2, equal to the right-hand side of equation (|6J), and 
Qr — = fic(l — 



further written as 



After some straightforward algebra, we find that this condition can be 

\LiJi 



1 



(47r)2-5 (if). 



[l+UJe 



(65) 



where ac denotes a critical value for the electrical conductivity, equal to the right-hand side of 
equation ([65|); so, numerical study of the latter condition requires calculation of the electrical 
conductivity a of the model under consideration. 

Electrical conductivity of compact stars has been studied by several investigators (see, e.g.. 
Flowers &: Itoh 1976| ; Itoh et al. 1983, hereafter I&83; Itoh & Kohyama 1993; Itoh, Hayashi 



fc Kohyama 1993| , hereafter IHK93; potekhin 19"96| ; [Potekhin fc Yakovlev 19961 ; [Potekhin et~ar 



19991 ; iPotekhin 1999| , hereafter P99).When ignoring the influence of magnetic field on electrical 



conductivity, we can proceed to the calculation of a by constructing a software package based, e.g., 
on I&83 (eqs. [2], [4]-[6], [8]-[9], and data from Table 1) and on IHK93 (eqs. [5], [7], [9], [28], [35], 
[39], [41], [43]-[44], [47], and data from Tables 1-3). For model parameters chosen as in § 6, typical 
values concerning one-component plasmas (e.g., He plasma with atomic number Z = 2 and mass 
number ^ = 4, C plasma with Z = 6 and A = 12, O plasma with Z = 8 and A = 16, etc.) are 



- 22 - 



a ^ 10^^ 10^" s ^, calculated at the base of the surface zone, and ~ 10^^ s ^ . So, the condition 
^) is easily satisfied by our model, since a joc^ 10^ 10^. 



The main difference between the above study for a perpendicular rotator and the study for a 
turning over model (made in the framework of MMM) is that in the first case differential rotation 
is removed exclusively when the model is already perpendicular rotator; while in the second case 
differential rotation is progressively removed as the turn-over angle increases, and when the model 
becomes perpendicular rotator its rotation is already rigid. So, there is need to examine in detail 
the reasons which lead to angular momentum mixing in both models. 

In particular, angular momentum mixing takes place mainly due to the action of the Alfven 
speed along the Wmxt direction (remember that changes most rapidly along exactly this di- 



rection), 'VAvj\mx(\ (generalization as in Chandrasekhar 1981 , Chapt. IV, eq. [60]), which carries 



angular momentum via hydromagnetic transverse Alfven waves propagating along the poloidal mag- 
netic component on this direction. In the perpendicular rotator, the poloidal magnetic component 
on the Wmxt direction, inducing VA-cu[mxt]j is the Zmxr-component Hpz[mxr]- Remark that the mag- 
netic field remains axisymmetric in the mixer's system of reference; thus, equations (^)-(^) for the 
magnetic field are expressed in terms of the mixer's coordinates, imxr = S,mxt(,= for simplicity), 
and 



^mxr \ 



{j^mxt, ^x) = ^rnxt l^x + V ^ ~ V ^ ~ " (^^) 



In the starting model, on the other hand, the poloidal magnetic component inducing V^ro[mxt] is 
the ^nixr'Component -f^pro[ma;r]- Generally, in any intermediate configuration (oblique rotator) the 
corresponding poloidal magnetic component, ^^pro[mxt]) is given by 

-^pm[mxi] (C) ^mxt ) l^x) ~ '^p ' ^xjj[mxt\ ~ 



Hpzu[mxr]{Ct l^mxriymxti ^x)) ~^ pz[mxr\i,^ i ^mxr{l^mxti ^x))\f^ ^ 

(67) 



,2 . 

X ' 



so, when i/^. (perpendicular rotator) we verify the first case, and when z^-^ ^ 1 (starting model 
with X ^ 0°) verify the second case. 

Now, the difficulty with the starting model has to do with the fact that, for model parameters 
chosen as in § 6, 



-^pz [mxr] 
-^p-uulmxr] 



10: 



in other words, angular momentum mixing in the starting model should be carried out by the 
"weak" poloidal magnetic component, -f^pro[mxr]- Then, after turn-over has started (assuming that 
Hpvj[mxr] can indeed drive angular momentum mixing), the magnetic field succeeds in progressively 
removing differential rotation for any intermediate configuration, since (1) angular momentum 
mixing is then carried out by the component -ffpro[mxt] (> -f^profmzr])) and (2) 0,3c{t) and Fr{t) 
decrease as t increases (thus, as seen by eq. o"c decreases as t increases). 
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On the other hand, magnetic fields in compact stars complicate electron transport making 
it particularly anisotropic (see, e.g., P99, § 1). As a result, electrical conductivity becomes a 
second-rank tensor (P99, eqs. [16]-[18]) with prevailing components the longitudinal electrical 
conductivity, a"||, acting along the resultant field lines, and the transverse electrical conductivity, 
(T_L, acting across the resultant field lines. To take into account such effects, we introduce a two-fold 
generalization. First, the electrical conductivities of interest are now written as 

^7,1 = (H) = c7|| (Ht + Up) ~ {Ut) , (69) 

CT± = ai_ (H) = ai_ {Ht + Up) ~ (Ht) (70) 

(magnetic fields as in eqs. since the ratio Ht/Hp is large over the whole star, except 

for two confined regions near the center and near the boundary; typical values for our model are 
Ht/Hp ~ 10^ -i- W^. Now, in the starting model, angular momentum mixing is mainly carried out 
by hydromagnetic Alfven waves propagating along the field lines of the poloidal magnetic field, 
that is, across the field lines of the resultant magnetic field, H, which, in turn, does almost coincide 
with the toroidal field, Hj. Remark that, for the starting model, the corresponding Alfven waves 
propagating along the resultant field lines do not contribute to the mixing process since, by this 
way, they are trapped to travel on cylinders of constant angular velocity. Accordingly, angular 
momentum mixing is mainly monitored by the transverse electrical conductivity of the toroidal 
magnetic field, a±{Ht). Thus, in summary, in the starting model angular momentum mixing is 
powered by -ffpro[m.2:t] (as discussed above) and its efficiency is monitored by cr±{Ht) (as discussed 
here) . 

Second, the quantities p2 and cTc become now functions of position and time (in the rest of 
this section, we shall use the mixture's system of reference; and, for simplicity, we shall denote the 
coordinate i^mxt by u alone). 



Hl_,^ ^^^{^,u,u^{t)) 
47r^(C, z.) nlit) iv^siC, u),Frit)) a2s2(e, u) ' 



' ^ {47r)^-^nsc{t)u:{s{^,iy),Fr{t))a^s^{^,u)pl-^{^,u,t)' ^ ^ 

where, in addition, we have substituted the poloidal field Hp by its component Hp^y^^t] along 
the 

'^mxt direction, since we are interested for the corresponding Alfven speed VA-nj[mxt]j which, as 
discussed above, has maximum efficiency in interchanging angular momentum. 



On the basis of such a two-fold generalization, the condition ( pq ) is now written as 

' Cr±iHt(C, Vraxr{v, V^{t))))\ 



(yc{i,l',t) / 



> 1 . (73) 



To compute the transverse electrical conductivity cr±{Ht), we use the Fortran Code CONDUCT, down- 



loaded from the site http://www.ioffe.rssi.ru/astro/conduct. This computer code has been 
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written by Potekhin and, actually, is the code that implements the formulae set up in P99 (see note 
in P99, last paragraph of § 6). Typical values for the starting model (with model parameters chosen 
as in § 6), when setting x(0) = Xini = vr/128 ~ 1.41° (which gives in turn z^;^(0) = = 0.9997) 

and when taking as typical case a carbon plasma, are {cr±{Ht)/ac) ~ 10'^. 

Such values, although apparently satisfying the condition ([73D, reveal a critical issue concern- 
ing the starting model. In particular, it seems that the toroidal magnetic field must be strong 
enough so as to induce an adequate dynamical asymmetry for the turn-over to take place; while, 
at the same time, it must be weak enough so as to permit interchanging of angular momentum 
via hydromagnetic Alfven waves propagating across its field lines. The difficulty is that, on the 
one hand, dynamical asymmetry increases with the toroidal magnetic field (i.e., with the magnetic 
perturbation parameter h) and, on the other hand, efficiency in interchanging angular momentum 
decreases as h increases (that is, (T±{Ht) decreases as h increases). In view of this remark and to 
the extent that our problem concerns mainly the behavior of the starting model, the right title for 
the present section seems to be: "is the toroidal magnetic magnetic field strong enough to induce 
dynamical asymmetry and weak enough to permit angular momentum mixing?" . 



8. Results and discussion 

In our computations, we first choose M = O.89M0 and = 1.00 (physical differential rota- 
tion). Taking then as typical case a carbon plasma with luminosity £ ~ lO^^ergs"^ (which gives 
in turn interior temperature 1^ ~ 3.97 x lO'^K, calculated by eq. |jl5|), and assuming for the 
starting model x(0) = Xini = vr/128 ~ 1.41° (which gives in turn ly^i^O) = v^lini] = 0.9997), we 
compute the quantities (/11//33 — 1) x 10, which is a measure of the dynamical asymmetry, and 
W±{Ht)/o'c) X 10^^, which is a measure of the degree that toroidal magnetic field permits angu- 
lar momentum mixing, in terms of h. The computations are performed with angular momentum 
Lxx - {L{Pnow) + ix)/2 (symbols as in §§ 3-4). 

Figure 1 shows these two quantities as functions of h. We remark that an "optimal value" 
for the magnetic perturbation parameter is /i ~ 0.08 in the sense that at this value the ratio 
{(T±{Ht)/ac) X 10~^ ~ 1.1 is still greater than unity, while at the same time the quantity (/11//33 — 
1) X 10 ~ 0.3 verifies that dynamical asymmetry has been already set up. We therefore take 
hDR = huB. = h = 0.08 for all the subsequent computations. 

For the constant-mass sequence [M = O.89M0, h = 0.08, Fr = 1.00], Figure 2 shows the graph 
of the function P = P{L). The point a{L = 0.979 x lO^^ergs, P = 52.091s) represents an almost 
spherical model due to counterbalancing of the effects of rotation and poloidal field (both tending 
to derive oblate configurations) with the effects of toroidal field (tending in turn to form prolate 
configurations); the effective ellipticity of this model is |ee//| = \{^av — Cp)\/Cav < W~^. The 
point b(L = 1.615 x 10^^ergs,P = 33.080 s) is a model with period equal to Pnow] so, L{Pnow) = 
Lf, = 1.615 X 10^^ ergs. The point c{L = 2.519 x 10^^ ergs, P = 22.150s) is the model for which 
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111 = (dynamical asymmetry holds left to this point); so, Lx = Lc = 2.519 x 10^^ ergs. The 
point d{L = 7.904 x IQ-^^ ergs, P = 12.634s) is the model of minimum period, Pmin = 12.634 s; so, 
L{Pmin) = = 7.904 X 10''^ ergs. This point defines two branches on the graph; in particular, on 
the branch left to d, secular evolution of the star induces a spin-down; while on the branch right 
to d (that is, evolution from e to d), secular evolution induces a spin-up. As in a recent paper 
( Geroyannis fc Papasotiriou 2000| , hereafter GPOO, § 4), we assume that secular evolution of the 
star is initiated at a point left to d and almost coinciding with it; eventually, this assumption means 
that the star does not undergo any secular spin-up during its evolution. The secular spin-down 
time rate, Pgec, is involved in the relation (GPOO, eq. [3]) 

/ 3c3 \ ^ 

Bs^i—IRjP,,,Psecj , (74) 

which yields as typical values for our model Psec ~ 10^^ -i- 10~^ Pnow with angular momentum lying 
in the interval [Lb,Lx]; in other words, the secular spin-down time rate is totally screened by the 
spin-down time rate observed now, as suggested in the first paragraph of § 3. Finally, the rightmost 
point e{L = 16.553 x lO^^ergs, P = 16.582 s) represents a critically rotating model with ratio of 
kinetic to gravitational energy r = T/|I^| ~ 0.14 and effective ellipticity Ce// — 0.39. 

Numerical solution of the system ( |60|) gives as optimal values for our model L^x = 2.154 x 
10^^ ergs, Bg = 2.525 x 10^ G, and 6 = 3.864. Then Figure 3 is a magnification of the graph 
between points b and c of Figure 2 together with the corresponding graph of the function Prr = 
Prr{L) (i.e., the period of the perpendicular rotator represented by the rigidly rotating terminal 
model). The arrow in this figure represents the "optimal nonaxisymmetric transition" from the 
axisymmetric differentially rotating starting model (aligned rotator) with Lxx = 2.154 x 10^^ ergs 
and Pxx = 24.949 s to the nonaxisymmetric rigidly rotating terminal model (perpendicular rotator) 
with Lrji = Lxx and Prr = 40.096 s. 

As discussed in § 3, Prr = 27r/r2i, where $7i is given in turn by equation (^ with the 
clarification that lRR^ii{h, L) = Ar{Ii, L)Iii{h) (note that, for the needs of the present section, we 
shall write explicitly the angular momentum L as second independent variable in this functional 
relation). The correction Ar{Ii,L), so-called rigid rotation amplification ratio, is calculated by 
the algorithm developed in § 5. Figure 4 shows the "noncorrected value" ( "noncorrected" in the 
sense that we attempt to approximate by this value the corresponding value of the perpendicular 
rotator) /ii(/i) = 1.323 x lO^'^gcm^ and the corresponding "corrected values" Ar{Ii, L)Iii{h) for 
several angular momentum values lying in the interval [Lh,Lx]. The arrow in this figure shows 
the corrected value for L = Lxx, ^^{Ii, Lxx)Iii{h) = 1.375 x lO^^gcm^; i.e., Ar{Ii, Lxx) = 1-039, 
which is equivalent to a percent correction ~ 4% relative to /ii(/i). Now, as discussed in § 4, 
the leading error involved in such calculations is due to the approximating relation (^). Figure 5 
gives lRR^ii{h, L) as calculated (1) by 2DCIT, and (2) by the approximating relation (pS]). This 
figure reveals a maximum error ~ 0.25% relative to the 2DCIT value, which means in turn that 
the leading error involved in the calculation of the correction Ar{Ii,L) is ~ 20 times less than the 
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(percent) correction itself; consequently, Ar{Ii,L) is a reliable estimate. 

Taking into account (1) the fundamental assumption of TOV scenario that the spin-down time 
rate remains constant during the turn-over and equal to the value observed now, 

P = constant = Pnow = 5.64 x 10"^^ ss"^ , (75) 

and (2) the above given optimal values (for L^x, Bg, and 5), we find t_D_RD = trov = 8.517 x 10^ yr, 
DEdrd = DEtov = 1.008x10^^ erg, tnaw = {Pnow- Pxx) / Pnow = 4.572xlo6yr, s.ud Drovitnow) = 
DEdrd /toRD = 3.754 x lO^^ergs"^; this value is ~ 30 times less than the classicaly estimated 
spin-down power (eq. Q) and thus the spin-down problem can be drastically simplified. 

Furthermore, assuming as above that secular evolution is initiated at the point d of Fig. 2 and 
proceeds left to it, the following estimate for the secular timescale, tsEC^ can be given 

P — P 

. ^ XX ^ mm i^n\ 

tsEC - ,p X , (76) 

\ ' [Lxx,Ld] 

where the subscript [Lxx , L^] denotes the interval of angular momenta over which the average 
of the secular spin-down time rate is calculated. In the latter calculation, the magnetic flux is 
assumed invariant, fxx = fdi so, when we know the flux at L^x, fxx = BsRav{Lxx), we can find the 
average surface poloidal field holding at d, Bd = BgRaviLxx) / Ravi^d) and then, by using equation 



(74), the corresponding secular spin-down time rate Pd- The same can be done for any angular 



momentum lying in the interval [Lxx,Ld]. We thus find (Psec) ~ 6.1 x 10 ^^ss ^ and 

tsEC — 6.4 X 10® yr; thus, the secular timescale tsEC results ~ 2 orders of magnitude greater than 
the turn-over timescale tTOV- 

Next, Figure 6 gives the graph of the function x = x(^) iii the time interval [0, tTov] where, 
obviously, the time t = Q corresponds to the beginning of the turn-over; the arrow head in this figure 
shows the "current value" Xnow = x{tnow) — 71.6°. Furthermore, Figure 7 gives the graph of the 
function Fy = Fr{t), with the arrow head pointing to the current value F^^now] = Fr{tnow) — 0.32. 
Likewise, Figure 8 gives the graph of the function Dtov = Drovit) and the arrow head shows the 
current value DTov\now] = DTov{tnow) - 3.8 x lO^^ergs"^. 



At this point, it is worth remarking that the time dependent scheme adopted in § 6 (eq. [g^) is 
an alternative way to (re)calculate some significant time dependent quantities. This can be done by 
assuming an "approximate rigid body description" of our model. In particular, the time dependent 
angular momentum components Li{t) and L^^t) along the principal axes Ii and Is, respectively. 



can be written as (see, e.g., Goldstein 1972, eq. [5.20] for the rigid rotation analog 



Li(t) = Lxx,Jl - = CL{t)hi{mic{t) {uj),{t) , (77) 
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where CL{t) is a correction due to differential rotation and is the average over the coordinate 

s of the differential rotation function uj at the time t. For t = we have ii(0) = and 



^3(0) = Lxx = CLiO)lDR,33ih,L,x)^3M (u;),(0) , 



(79) 



where f^3c(i) = Q{0,{Lxx),t) {Q is the symbol for the time dependent scheme [p2|) and thus 
Q3c(0) = Q{Lxx); also, (w)^(t) = {uj{Fr{t)))^ = (uj{Q{F^[ini],t)'j ^ and thus (u;)^(0) = {^uj{F,.[ini]] 
So, the correction Cl(0) is given by 



Cl(0) 



lDB,33{h, Lxx)^{Lxx) \^{Fr\ini]) 



(80) 



Furthermore, CL{tTov) = 1 since at the time troy rigid rotation has been established throughout 
the star. For the intermediate configurations, we assume the time dependent scheme (^2|), 



CL{t)c^l + Q[[CLiO)-lit 
Thus, for the angular velocity components we can write 

Lxx^Jl- v^{tf 



(81) 



Olc(t) 
^3c{t) 



CL{t)Iu{t) Hsit) 

Fxxi^x(^) 



(82) 



(83) 



Consequently, the resultant angular velocity along the z axis (keep in mind, however, that the 
resultant angular velocity inclines at a small angle 7 with the z axis), 0,zcit), is given by 



(84) 



The above involved time dependent moments of inertia /ii(t) and /33(t) are approximated by the 
relations 

hi it) ^ Ar ih ,Lxx)In{h)-Q( [Ar (h , Lxx)hi (h) - Idrm {h,Lxx)lt^ , (85) 

(86) 

For t = we get /ii(0) = lDR,ii{h,Lxx) and 133(0) = lDR,33ih, Lxx); while for t = trov we find 
hiitrov) = Ar{h,Lxx)Iii{h) and l33{tTOv) = l33{h). 



I33it) ^ l33{h) + Q [lDR,33{h,Lxx) " h3ih)],t 



For the rotational kinetic energy T(t) we can write (see, e.g., Goldstein 1972, eq. [5.21] for the 
rigid rotation analog) 



in{t){nut) (^),(t) ]\i33{t)(nut) (^>,(t)' ' 



(87) 
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where Cxit) is a correction due to differential rotation. For t = we have 
T(0) =r,.,. = ^CT{0)lDRMh,L^.)(^^3c{0) Ms(o))' 



which gives 

Ct(0) = _21^2_ ^gg^ 



lDR,33{h, Lxx) (^^{Lxx) (^to{Fr[ini]) 

Furthermore, CxitTOv) = 1 since at the time tTOV rigid rotation has been estabhshed throughout 
the star. For the intermediate configurations, we adopt the time dependent scheme (52), 

CT{t)c:^l + Q(^[CT{0)-l],ty (90) 

Now, we are ready to (re)calculate some significant quantities on the basis of the approximate 
relations (|77|)"-(|90|). In particular, the central period P{t) is written as 

where ^zc{t) is given by equation (U). On the other hand, however, the fundamental assumption 
(|75r) leads to the linear relation 



Pit) = Prb - (Prr - Pxx)Mt) , (92) 
where ^'^(t) is defined by equation (p^). To distinguish between the values (pT])-(92) of P{t), we add 



the subscript FB to the first one, PpB = '^^^/^zc^ to give emphasis on that this is a "feed-back" value 
for P{t). In order for the model to retain its inherent consistency, these two values should be close 
to each other, (P(t) - PpB(t)) x W'^/P{t) < tol, where tol is a reasonable percent tolerance. Figure 
9 shows (1) the function P{t), which results from the assumption ([f5|), and (2) the function PpBit) 
of which the time derivative Pps is a feed-back value for the spin-down time rate, slightly deviating 
from the initial assumption ([75|). This figure reveals a maximum percent difference ~ 2% relative 
to P{t), which, however, when taking into account the approximate character of equations (|77|)- 
(pO|)), seems to be insignificant. Furthermore, the corresponding average (^PpB^^ — 5.6 x lO"^"^ ss~^ 

over the whole turn-over time deviates less than 1% from Pnow In addition. Figure 10 shows 
(1) the central angular velocity ^zc{t) calculated with constant P; (2) its counterpart ^zc[fb] 
calculated with feed-back P\ (3) the central angular velocity ^^dt) calculated with constant P\ 
(4) its counterpart ^^c[fb] calculated with feed-back P; (5) the central angular velocity ^ic{t) 
calculated with feed-back P; and (6) the angular velocity Q-z{t) of the progressively established rigid 
rotation in the mixture configuration approximated by the relation Q-z{t) ~ {ytzc{t) — ^'ic{t)) {^)sit)- 

Next, Figure 11 shows the rotational kinetic energy T(t) of the model, calculated (1) by the 
relation T{t) = Txx — Jo Drovit) dt (where DTov{t) is given by eq. |59| with ixov substituted 



by t), and (2) by the relation (^) which, as discussed above, represents the corresponding feed- 
back value TFB{t). As seen in this figure, the extremum percent difference is ~ —1.8% relative 



- 29 - 



to T{t) near t = A x 10^ yr. Since, however, in our model the quantity T is of higher significance 
than T itself, it is interesting to remark that these two graphs exhibit almost equal time rates; 
in addition, their average time rates over the whole turn-over time are almost equal: (Tfb 

f) = -3.754 X 10^2 erg s-^ 



t ~ 



t 

Figure 12 shows the angle 7(t) between the instantaneous angular velocity axis and the in- 
variant angular momentum axis. As it can be verified after some elementary algebra, this angle is 
given by 

Mt) = ^ + ^] (93) 

The arrow head in this figure shows the current value 'jnow = jitnow) — 3100 arcseconds. The 
maximum value of 7 is ~ 3500 arcseconds near t = 3.5 x 10^ yr; thus, as emphatically suggested in 
the foregoing sections, the angle 7 remains less than 1° during the whole turn-over. 

Finally, Table 1 gives a summary of calculations for both the cases F^Uni] = 1.00 and 0.50. 
. ... . . .... ^ . .0 ... . ........ . 

the spin-down problem in AE Aquarii. Concluding, we find interesting to emphasize once again on 
the remark made in the Abstract: "an observed large spin-down time rate does not always imply 
a large spin-down power" . 



The TOV scenario has resulted as a product (contributed by the author) within the frame- 
work of a research project supported by the Research Commitee of the University of Patras (C. 
Caratheodory's Research Project 1998/1932). It is a pleasure for the author to acknowledge his col- 
leagues in this project Dr. Florentia Valvi and Dr. Panagiotis Papasotiriou for helpful discussions 
and comments on the manuscript. 
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Table 1. Summary of calculations regarding TOV scenario' 



Parameter'^ 



1.00 



0.50 



Optimal Lxx 
Optimal Bs 

Optimal S (dimensionless) 
Magnetic flux, / 

Rigid rotation amplification ratio, Ar{Ii, L^x) (dimensionless) 
Secular timescale, tg ec (yr) 
TOV & DRD timescales, txov = toRD 
TOV & DRD timescales in units of P^x 
Present TOV time, tno™ (yr) 
Present spin-down power, DTOv(tnow) 

t(ofcq. 1^ ) 



(yr) 

(number of revolutions) 



inRD 



Ratio 



Power dissipated due to electron and ion viscosity /ie+i> 
Present turn-over angle, Xnow (arcdegrees) 
Present reduction factor, iV[noiD] (dimensionless) 
Present central angular velocity about the 1^ axis, Qsdnow] 
Present angular velocity - angular momentum angle, -f„ow (arcseconds) 



2.154(-|-49)'' 

2.525(-|-06) 

3.864(-|-00) 

1.105(-l-03) 

1.149(+24) 

1.039(+00) 

6.392(-|-08) 

8.517(-|-06) 

1.077(-l-13) 

4.572(-|-06) 

3.754(4-32) 

2.664(-|-01) 

3.467(-|-30) 
7.160(-|-01) 
3.156(-01) 
7.849(-02) 
3.102(-|-03) 



2.272(-|-49) 
1.018(-l-06) 
7.160(-l-00) 
1.001(-l-03) 
4.663(-l-23) 
1.043(-l-00) 
7.621(-|-09) 
4.873(-|-06) 
5.208(-l-12) 
2.009(-l-06) 
2.218(-(-32) 

4.509(-l-01) 

1.201(-l-30) 
7.071(-l-01) 
1.652(-01) 
7.036(-02) 
3.294(-|-03) 



Calculations for the axisymmetric differentially rotating starting model (aligned rotator) 



Central period, Pxx 
Rotational kinetic energy, Txx 

Moment of inertia along the Ii axis, I DR,i\{h, Lxx) 
Moment of inertia along the /s axis, lDR,33{h, Lxx) 
Average radius, Rav 

Radius (spherical approximation), Rsph 
Volume (spherical approximation), Vgph 
Average density (spherical approximation), (^Qsph^ 
Average pressure (spherical approximation), (^Psph') 
Average sound speed (spherical approximation), (^a^p^^ 
Moments of inertia (spherical approximation), IiHsph] = 
Dimensionless radius (spherical approximation), 
Dimensionless base of the surface zone (spherical approximation) , 
Average surface poloidal magnetic field, Bs 
Proper poloidal magnetic parameter, /9j 
Average surface toroidal magnetic field, {Ht{^B, >^)) ^ 
Surface Alfven speed (spherical approximation), 



'22[sph] — ^33[sph] 



As[sph] 

Surface Alfven time (spherical approximation), tyis[sph] = 

^A[sph] 



Average Alfven speed (spherical approximation), Vj 



Average Alfven time (spherical approximation), tyi[spfi] 
Maximum poloidal magnetic field, Hp^^^^^^ 
Average poloidal magnetic field, {Hp) 

Average z-component of the poloidal magnetic field, {Hpz) 
Average ro-component of the poloidal magnetic field, {Hp^) 
Maximum toroidal magnetic field, -ffjf^^j.] 
Average toroidal magnetic field, {Ht) 
Ratio f2c/r2e = l/oJe 



^sph/^As[sph] 
_ <Hp> 

Rsph/VA[sph] 



2.495(-l-01) 
1.789(-l-48) 
1.387(-l-50) 
1.376(-l-50) 
6.744(-|-08) 
6.525(-l-08) 
1.164(-|-27) 
1.490(-l-06) 
4.450(4-23) 
5.466(4-08) 
1. 285(4-50) 
3.711(4-00) 
3.628(4-00) 
2.525(4-06) 
2.011(4-04) 
2.810(4-09) 
6.501(4-03) 
1.004(4-05) 
8.618(4-03) 

7.571(4-04) 
1.172(4-08) 
3.729(4-07) 
3.461(4-07) 
7.271(4-06) 
7.760(4-11) 
3.198(4-11) 
3.625(4-00) 



2.951(4-01) 
1.904(4-48) 
1.391(4-50) 
1.384(4-50) 
6.767(4-08) 
6.527(4-08) 
1.165(4-27) 
1.487(4-06) 
4.439(4-23) 
5.463(4-08) 
1.286(4-50) 
3.710(4-00) 
3.627(4-00) 
1.018(4-06) 
4.985(4-04) 
2.811(4-09) 
2.621(4-03) 
2.490(4-05) 
3.475(4-03) 

1.879(4-05) 
4.722(4-07) 
1.502(4-07) 
1.394(4-07) 
2.929(4-06) 
7.750(4-11) 
3.194(4-11) 
2.125(4-00) 



Calculations for the nonaxisymmetric rigidly rotating terminal model (perpendicular rotator) 



Period, Pru 

Rotational kinetic energy, Tnn 



4.010(4-01) 3.817(4-01) 
1.688(4-48) 1.870(4-48) 
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Table 1 — Continued 



Parameter'' 1.00 0.50 



Moment of inertia along the /i axis, Ijifi Lxx^ 


1.375 


+50) 


1.380(+50) 


Moment of inertia along the I3 axis, 133(h) 


1.283 


+50) 


1.283(+50) 


Average radius, Rav 


6.574 


+08) 


6.574(+08) 


Radius (spherical approximation), Rsph 


6.409 


+08) 


6.409(+08) 


Volume (spherical approximation), Vgpfi 


1.103 


+27) 


1.103(+27) 


Average density (spherical approximation), (^Qsph') 


1.600 


+06) 


1.600(+06) 


Average pressure (spherical approximation), {^Psph') 


5.012 


+23) 


5.012(+23) 


Average sound speed (spherical approximation), ^a^pft^ 

Moments of inertia (spherical approximation), Iii[sph] = ^22[sph] = ^33[sph] 


5.598 


+08) 


5.598(+08) 


1.255 


+50) 


1.255(+50) 


Dimensionless radius (spherical approximation), 


3.730 


+00) 


3.730(+00) 


Dimensionless base of the surface zone (spherical approximation), 


3.650 


+00) 


3.650(+00) 


Average surface poloidal magnetic field, Bs 


2.658 


+06) 


1.079(+06) 


Proper poloidal magnetic parameter, /3J 


1.935 


+04) 


4.765(+04) 


Average surface toroidal magnetic field, {Ht{is,^)) ^ 


2.762 


+09) 


2.762(+09) 


Surface Alfven speed (spherical approximation), V^s[3p/i] 


6.842 


+03) 


2.778(+03) 


Surface Alfven time (spherical approximation), 4As[sph] 


9.368 


+04) 


2.307(+05) 


Average Alfven speed (spherical approximation), V^[sph] 


9.146 


+03) 


3.714(+03) 


Average Alfven time (spherical approximation), tA[sph] 


7.007 


+04) 


1.726(+05) 


Maximum poloidal magnetic field, -ffp[maa;] 


1.299 


+08) 


5.274(+07) 


Average poloidal magnetic field, {Hp) 


4.101 


+07) 


1.665(+07) 


Average z-component of the poloidal magnetic field, (Hp^) 


3.806 


+07) 


1.545(+07) 


Average ro-component of the poloidal magnetic field, (Hp^) 


8.004 


+06) 


3.250(+06) 


Maximum toroidal magnetic field, H^^^^j 


8.252 


+11) 


8.252(+ll) 


Average toroidal magnetic field, (Ht) 


3.387 


+11) 


3.387(+ll) 


Ratio Cc/f^e = 1/ti^e 


1.000 


+00) 


1.000(+00) 



''Basic model: M = 0.89 M©; h = 0.08; k = cos 45° = 0.7071; Xini = 7r/128 ~ 1.41°; ly^ii^i] = 0.9997 
''In cgs units, unless stated otherwise 

Parenthesized numbers denote powers of 10 
''Typical values for a carbon plasma with £ = 10^^ ergs"', Ts ~ 3.97 X 10^ K, Qs — 1-2 X 10* gcm~^ 
"Involved parameters: Ap = 1; A^„^ = 10; A.^nv = 2 X 10 (for details, see § 6) 



'Details as in 



Durisen 1973 



TOV ' TOV 

3, especially eqs. (4), (6) 
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Fig. 1. — The quantities {(^±{Ht)/ac) x 10~^ and (/11//33 — 1) x 10 as functions of h. The double- 
head arrow has base at h = 0.08 and shows for the former and latter quantities their corresponding 
values 1.105 and 0.313, respectively. This is an "optimal value" for the magnetic perturbation 
parameter h in the sense that at this value the former quantity is still greater than unity, while at 
the same time the latter quantity shows that dynamical asymmetry has been already established. 

Fig. 2.— P vs. L for the constant-mass sequence [M = 0.89MQ,h = 0.08, i^^ = 1-00]. Point 
a{L = 0.979 X 10^^ ergs, P = 52.091 s): almost spherical model due to counterbalancing of the effects 
of rotation and poloidal field with the effects of toroidal field. Point b{L = 1.615 x 10^^ergs,P = 
33.080s): model with period equal to Pnow', so, L{Pnow) = Lh = 1.615 x 10'*^ ergs. Point c(L = 
2.519 X lO'^^ergSjP = 22.150s): model for which In = I33 (dynamical asymmetry holds left to 
this point); so, 1^ = 1^ = 2.519 x lO'^^ergs. Point d{L = 7.904 x lO^^ergs,^ = 12.634s): 
model of minimum period, Pmin = 12.634s; so, L{Pmin) = Ld = 7.904 x 10^^ ergs. Point e{L = 
16.553 X 10^^ergs,P = 16.582 s): critically rotating model with ratio of kinetic to gravitational 
energy r = T/\W\ = 0.140 and effective ellipticity ee// = 0.388. 

Fig. 3. — Magnification of the graph between points b and c of Fig. 2 together with the cor- 
responding graph of the function Prr = Prr{L), i.e., the period of the perpendicular rotator 
represented by the rigidly rotating terminal model. The arrow shows an "optimal nonaxisymmet- 
ric transition" from the axisymmetric differentially rotating starting model (aligned rotator) with 
Lxx = 2.154 X 10^^ ergs and Pxx = 24.949s to the nonaxisymmetric rigidly rotating terminal model 
(perpendicular rotator) with Lrr = L^x and Prr = 40.096 s. 

Fig. 4. — The "noncorrected value" Iii{h) = 1.323 x lO^'^gcm^ ( "noncorrected" in the sense that 
we attempt to approximate by this value the corresponding value of the perpendicular rotator) and 
the corresponding "corrected values" Ar{Ii, L)Iii{h) vs. L in the interval [Lf,, Lx]. The arrow shows 
the corrected value for L = Lxx, Ar{Ii, Lxx)hi{h) = 1.375 x 10^'^gcm^; i.e., Ar{Ii, Lxx) = 1-039, 
which is equivalent to a percent correction ~ 4% relative to /ii(/i). 

Fig. 5. — lRR,ii{h, L) calculated (1) by 2DCIT and (2) by the approximating relation ( p5| ) vs. L in 
the interval [Lb, Lx\. This figure reveals a maximum difference ~ 0.25% of the latter relative to the 
former values, which means in turn that the leading error involved in the calculation of At.{Ii,L) 
is ~ 20 times less than the (percent) correction itself. 

Fig. 6. — The turn-over angle x vs. t in the interval [0, troy]- The arrow head shows the "current 
value" Xnow = x{tnow) = 71.604°. 

Fig. 7. — The differential rotation reduction factor Fy vs. t in the interval [0,tToy]- The arrow 
head points to the current value -^^[noiu] = Fr{tnow) = 0.316. 

Fig. 8. — The power Dtov dissipated due to turn-over vs. t in the interval [0,tToy]- The arrow 
head shows the current value Drpovinow] = DTovitnow) = 3.754 x lO^^ergs"^. 
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Fig. 9. — The central period P{t) as it results (1) from the assumption ([7^) and (2) from the 
function PpBit) (eq. of which the time derivative PpB is a feed-back value for the spin-down 
time rate, slightly deviating from the initial assumption ([75[). This figure reveals a maximum 
percent difference ~ 2% of the latter relative to the former values; however, when taking into 
account the approximate character of equations ([7^)-(|90[), this difference seems to be insignificant. 

Fig. 10. — (1) the central angular velocity i^zc{t) calculated with constant P; (2) its counterpart 
^zc[FB] calculated with feed-back P; (3) the central angular velocity i^sdt) calculated with constant 
P; (4) its counterpart f^sc^FB] calculated with feed-back P; (5) the central angular velocity 0,ic{t) 
calculated with feed-back P; and (6) the angular velocity i^zit) of the progressively established rigid 
rotation in the mixture configuration approximated by the relation 0,z{t) ~ (i}zc{t) — ^3c{t)) {^)si^)- 

Fig. 11. — The rotational kinetic energy T calculated (1) by the relation T{t) = Txx — Jq DTOv{t) dt 
(where DTovit) is given by eq. |^ with Itov substituted by t), and (2) by the relation ( p7| ) which 
represents the corresponding feed-back value Tpsit). The extremum percent difference of the latter 
relative to the former values is ~ —1.8% near t = 4 x 10^ yr. However, in our model the quantity 
T is of higher significance than T itself; it is therefore interesting to remark that these two graphs 
exhibit almost equal time rates; in addition, their average time rates over the whole turn-over time 
are almost equal: (Tfb^^ — ^ -3.754 x lO^^ergs"^ 

Fig. 12. — The angle 7(t) between the instantaneous angular velocity axis and the invariant angular 
momentum axis vs. t in the interval [0, tTOv]- The arrow head shows the current value ^now = 
litnow) = 3102 arcseconds. 
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